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1. Introduction 

The mass splitting in the meson and baryon system is an interesting topic in hadron spec- 
troscopy. It is related with the quark masses which are the fundamental parameters of the Standard 
Model. The mass difference of the proton and neutron is also crucial to the phenomenological 
model of nuclei, because it plays an important role in the neutron j8 decay process. These mass 
differences have already been measured by experimentalists to good accuracy. So they can be used 
as input to help theorists to determine the parameters in their phenomenological models or check 
the validity of the theory. 

The mass of the hadron is determined by both QED and QCD dynamics. For the QCD interac- 
tion, since the coupling constant is too large at the low energy regime, the perturbation theory is not 
applicable. People have developed all kinds of effective theories of QCD from the 1950s. They can 
be used to describe the hadron mass spectrum. For instance, the mass of the pseudoscalar meson, 
especially for the pion, is much smaller than the mass of the baryon. It can be explained by the 
spontaneous breakdown of chiral symmetry. When QED is included, the meson masses are split, 
but m 2 n+ — m 2 n Q ~ m 2 K+ — m 2 ^, which is explained by Dashen's theorem [jl|]. Partially quenched chi- 
ral perturbation theory(PQ^PT) gives a more general description of QCD phenomena: valence and 
sea quark masses are varied independently. It is important since the cost to generate the dynamical 
QCD configurations for a large lattice is still pretty high today. 

Recent developments in PQ^PT [§] allows us to incorporate the QED interaction too, so we 
can investigate the mass splitting which is related to the QED interaction. Since the QED interac- 
tion is much weaker than the QCD interaction, it will only provide a small shift in the hadron mass 
from its QCD value. But this shift is important to the splitting in the pion and kaon multiplet. 

Because the lattice calculation can give the hadron mass spectrum from first principles, it is 
used to test whether the QCD+QED interaction can actually reproduce the hadron spectrum. The 
pioneering work on this topic has been done in ^ with the quenched QCD and non-compact QED 
configuration. The development of Lattice QCD enables us to investigate the hadron spectrum 
with 2+1 flavor dynamical QCD configurations, which means we can include the effect of the 
degenerate up and down quarks and also the heavier strange quark in the sea quark sector. So our 
simulation result can be used to fit the PQ^PT formula. We can extract the unknown low energy 
constants(LEC's) from the fitting. After we know the LEC's, we can use the pseudoscalar meson 
mass as input to determine the physical value of the u, d, and s quark masses. The MILC group 
also reports their analysis of the mass splitting using staggered fermion ensembles in Basak's Talk. 

The lattice simulation is also used for the mass splitting in the nucleon system [Q]. We present 
preliminary results for the mass difference between the proton and neutron due to QED. 

2. Theoretical Background 

In chiral perturbation theory of QCD, the mass square of the pseudoscalar meson is determined 
by the spontaneous breakdown of chiral symmetry from SU (3)l ® SU (3)r®U (l)y to SU (3)y <g> 
t/(l)y. The Goldstone boson generated by this mechanism is the pseudoscalar meson octet. All 
of them are massless in the chiral limit. If the QED interaction is also considered, the symmetry 
group will be broken further to SU (2)y ®U{l) em ®U (l)y. The only Goldstone bosons left are the 
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neutral pion and kaons (the neutral pion can be considered a Golds tone boson if G{a 2 ~) terms are 
neglected[@]). 

Since we calculate with several valence quark mass and charge combinations, we use PQ^PT. 
Suppose we have a "kaon" which is composed by u and s quarks. Following the notation of 
index 1 denotes the u quark and 3 the s quark. The square of the meson mass is[^|] 



m 
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Xl3 + -=2- Wl 
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Yj are independent linear combinations of the LEC's K Er as defined in Eq. (48) of [I2J, and Xi 



is the LO meson mass, %i = 2B m h % i} = (Xi + Xj)/ 2 > ?13 
functions A,B,B\,R, and others can be found in [0]. 



q\ — q-i = q u — q s . The form of the 





QCD LEC's 


QED LEC's 


LO 


Bq,Fq 


C 


NLO 


jr jr jr jr 


C,Y\,Y 2 ,Y 3 ,Y 4 ,Y 5 



Table 1: LEC's in PQjPT 

Table |^ shows all of the LEC's which are necessary to determine the pseudoscalar meson mass 
at NLO. The QCD LEC's have been computed in[|7|], so the current simulation is used to extract 
the QED LEC's. Note that the definition of the decay constant used in ^ (followed here) is \fl 
smaller than the one in [^]. For the splitting in baryon system, we focus on the mass differences 
between the proton and neutron. Here, we study only the degenerate case, m u = m^. Then the 
leading contribution to the mass difference is proportional to a. 



3. Lattice Simulation 

Following [||], the lattice calculation employs combined QCD+QED gauge configurations. For 
the QCD configurations, we use the Nf = 2 + 1 flavor QCD configurations generated with the 
domain wall fermion and Iwasaki gauge actions by the RBC and UKQCD collaborations []|, |(| 0]. 
For the 16 3 x 32 lattice, we use light quark mass 0.01, 0.02 and 0.03 ensembles. For 24 3 x 64 
lattice, we use the 0.005, 0.01, 0.02 and 0.03 ensembles. The strange quark mass is fixed to 0.04 in 
all cases. Table || gives the details of the lattice ensembles used in this work. 

The QED configurations are generated in a quenched, non-compact form[]8|, 0]. Here we em- 
ploy the Feynman gauge instead of the Coulomb gauge which was used previously^]. One ad- 
vantage of the non-compact QED formalism is that the U (1) gauge potential can be chosen 
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lat 


^sea 


m va i 


Trajectories 


A 


meas 


I src 


16 3 


0.01 


0.01-0.03 


500-4000 


20 


176 





16 3 


0.02 


0.01-0.03 


500-4000 


20 


176 





16 3 


0.02 


0.01-0.03 


500-4000 


20 


176 





24 3 


0.005 


0.005-0.03 


900-8060 


40 


180 





24 3 


0.01 


0.01-0.03 


1460-5040 


20 


180 





24 3 


0.02 


0.02 


1800-3580 


20 


180 


0,32 


24 3 


0.03 


0.03 


1800-3580 


20 


180 


0,32 



Table 2: QCD gauge configuration ensembles generated by the RBC and UKQCD collaborations [|[ ^ ffa. 
A is the separation of the trajectory number in molecular dynamics time units. )3qcd = 2.13. The inverse 
lattice scale is a~ l = 1 .729(28 )GeV for both of 16 3 and 24 3 lattice. 



randomly with the correct distribution in momentum space, then Fourier transformed to coordinate 
space, so there are no autocorrelations in the ensemble. Another advantage is that there is no (lat- 
tice artifact) photon self-interaction in the action. To couple to the fermions, the non-compact 
field is exponentiated in the usual way, so the combined gauge field appears as an SU (3) matrix 
times a phase, U® CD (x)xU® ED (x). 

Since the QED interaction does not have confinement, it is possible that the finite volume 
can induce a significant systematic error. So, we do our simulation on both 16 3 and 24 3 lattice 
configurations with the same lattice spacing. Then we can investigate the finite volume effect in 
the mass spectrum. 



4. Numerical Results 

We can switch on and off the QED interaction by setting the charge e ^ and e = 0. We first 
extract the masses from wall source, point sink, correlation functions. We define the mass differ- 
ence, 8m 2 = m 2 (e ^ 0) — m (e = 0), and then fit 8m 2 with the PQ^PT formula, using the standard 
jackknife method to estimate the statistical error. 8m is averaged over +e and — e ensembles to 
decrease the noise in the signal [§]. We take the QCD LEC's from the St/ (3) PQ^PT fit in [§] since 
we use the same ensembles. We also use these for the 16 3 case because these two ensembles have 
the same lattice scale. Then we fit our data to obtain all of the QED LEC's, except Y\ because the 
sea quarks are not coupled to photons. Thus, Y\ = in this work. 

PQ^PT is based on chiral symmetry, which is not exact in the domain wall fermion formalism 
with finite L s . To lowest order, the violation of chiral symmetry can be counted as a small shift 
in the input quark mass, so we substitute m, + m res for all quark masses, where the residual mass 
m res is determined from the method in reference [Jsj] . Also following [^J, we add additional LEC's, 
C\ and C2, proportional to {q\ — q^) 2 and (q\ +qj) 2 , respectively, to account for chiral symmetry 
breaking due to finite L s at order &{a). These LEC's vanish when L s — ► 00. C\ mixes with the LEC 



C in Eq. (2.2), so it must be subtracted to obtain the physical LEC. C\ is measured from m res with 
e 7^ as described in [Q]. For the preliminary results presented here, this subtraction has not been 
performed. 
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By fitting the results for 5m 2 , we can get all of the QED LEC's, including the pure lattice 
artifact C2. The fit results are shown in Figure [j] along with the unitary data points (though all of 
the partially quenched data was used in the fit) and tabulated in Table |3| We perform fits with and 
without chiral logs. # 2 /dof is adequate in either case, except for 24 3 when all the quark mass data is 
fitted. If the heaviest sea quark points are omitted, x 2 becomes reasonable. The main effect of the 
logs is to significantly reduce the value of the charged meson splitting in the chiral limit. Note, the 
log terms vanish for the neutral mesons, and the neutral splittings should also vanish in the chiral 
limit (see [Jsj] for a detailed discussion). The non-vanishing chiral limit value in the plot is due of 
the C2 term, or the explicit chiral symmetry breaking induced by finite L s . 

16 3 unitary point mass square splitting 24 3 unitary point mass square splitting 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

m 2 m 2 
m ps m ps 

Figure 1: 16 3 (left panel) and 24 3 (right panel) lattice data points and fit results for the meson mass splitting. 
Lines correspond to fits with and without chiral logs. The data points in the plot correspond to ud, uu and dd 
mesons, respectively, from top to bottom. 





16 3 


24 3 


16 3 


24 3 


24 3 


Ref. \§ 


with log 


No 


No 


Yes 


Yes 


Yes 


Yes 


fit range 


0.01-0.03 


0.005-0.03 


0.01-0.03 


0.005-0.03 


0.005-0.02 


N/A 


Cx 10 6 


3.48(40) 


3.02(14) 


0.85(16) 


0.769(64) 


0.96(10) 


2.71 


Y 2 x 10 2 


1.47(3) 


1.51(1) 


1.56(3) 


1.54(1) 


1.51(2) 


0.53 


Y 3 x 10 3 


7.64(72) 


7.61(39) 


-1.08(45) 


-1.17(35) 


-2.79(43) 


2.65 


Y 4 x 10 3 


8.00(78) 


7.55(50) 


9.71(92) 


8.22(51) 


8.85(80) 


5.7 


Y 5 x 10 3 


1.48(23) 


2.43(12) 


1.05(24) 


2.23(12) 


2.35(19) 





C 2 x 10 3 


2.69(10) 


2.39(5) 


2.69(10) 


2.39(5) 


2.48(6) 


N/A 


X 2 /dof 


0.7(6) 


0.6(2) 


1.0(6) 


2.8(5) 


1.2(6) 


N/A 



Table 3: The QED LEC's from fits of 5m 2 . The pure QCD LEC's were taken from the SU(3) fit in [|]. The 
last column gives phenomenological results [Q]. 

From Table || we can see that the LEC's change slightly from 16 3 to 24 3 , but C and F3 change 
significantly when the the log terms are included. The effect on the other LEC's is not so drastic. 
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Also shown in Table |3| are the values from phenomenology which are roughly consistent with 
ours within an order of magnitude, though we emphasize that our values do not include systematic 
corrections due to finite volume, explicit chiral symmetry breaking, and non-zero lattice spacing, 
to name just three, and so should be taken as preliminary. 



Finite Volume Effect 
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Figure 2: Finite volume effect on the ud meson. Figure 3: The proton-neutron mass difference. 

Figure ^| shows the comparison of the charged meson splitting from 16 3 and 24 3 lattices, and 
the relative change is given in Table ||. The difference actually decreases at the smallest quark mass, 
0.01, which is why the finite volume effect shown in Figure ^ appears to go away in the chiral limit 
(linear fit). Interestingly, from Figure |TJ one can see that the neutral splitting is mostly insensitive 
to the finite volume, at least for the volumes studied here. We are still investigating these effects, 
and plan a detailed comparison with finite volume chiral perturbation theory [Q]. 



nival 


5m 2 (24 3 )-5m 2 (16 J ) 


Z Sm 2 (24 i )+Sm 2 {16 1 ) 


0.01 


0.1043(9) 


0.02 


0.1349(9) 


0.03 


0.1286(9) 



Table 4: Finite volume effect on the ud meson. 

Figure [| is the result of the mass difference between the proton and neutron when the QED 
interaction is included. If there is no QED interaction and m u = m^, then m n = m p , which is the 
result of isospin symmetry. Figure |3| shows the mass difference for the unitary points. We can see 
the proton is heavier than the neutron and the mass difference decreases with quark mass. The 
24 3 result is larger than the 16 3 result, once again signaling finite volume corrections. The baryon 
mass difference is noisier than in the pseudoscalar case. We plan to increase the statistics for the 
baryon mass spectrum to improve the precision on this result, and to compute the mass splitting 
coming from non-degenerate u and d quark masses, which is expected to switch the sign of the 
mass difference, in accord with Nature. 
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5. Summary and outlook 

Using 2+1 flavor QCD and non-compact QED, we have studied the electromagnetic mass 
splitting in hadrons. We have done our simulation on different sea quark mass ensembles and have 
used both 16 3 and 24 3 lattices to investigate the finite volume effect on the mass splitting. 

Preliminary fits appear to be consistent with NLO PQ^PT, including photons, for the param- 
eters in our study, except at the heaviest quark mass, 0.03, which is close to the physical strange 
quark mass. A more complete study of systematic errors, including finite volume and explicit chiral 
symmetry breaking effects, is in progress. Next we will use the LEC's obtained from our analysis 
to determine the physical, non-degenerate, u, d, and s quark masses. 
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